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ABSTRACT 


Statistical  procedures  for  shock  models  and  wear  processes  are 
considered  in  this  paper.  We  show  that  independent  gamma- 
Dirichlet  priors  are  conjugate  priors  when  sampling  from  these 
shock  models.  Bayes  rules  given  the  observations  are  computed. 
In  particular,  we  calculate  the  Bayes  estimates  of  the  survival 
probabilities  for  these  models.  We  show  consistency  of  the 
posterior  distribution  as  well  as  weak  convergence  of  the 
centered  and  suitably  rescaled  posterior  processes. 


BAYESIAN  NONPARAMETRIC  STATISTICAL  INFERENCE 
FOR  SHOCK  MODELS  AND  WEAR  PROCESSES 


by 

Albert  Y.  Lo 


1.  INTRODUCTION  AND  SUMMARY 


Let  a  device  be  subject  to  shocks  occurring  randomly  In  time  accord¬ 
ing  to  a  homogenous  Poisson  point  process  N  ■  (N(t)  ;  t  e  R}  with 
Intensity  X  .  The  1th  shock  causes  a  random  amount  of  damage  , 

1*1,  . . . ,  n,  ...  .  The  random  damages  are  assumed  to  be  Independent 
and  Identically  distributed  as  F  ,  where  F  Is  a  distribution  function 
supported  by  [(),<■)  .  (This  Implies  damages  are  never  negative,  also 
note  that  the  results  In  Section  2  do  not  depend  on  the  support  of  F  .) 

This  shock  model  have  been  studied  by  Esary,  Marshall  and  Proschan  (1973), 
Barlow  and  Proschan  (1975),  among  others.  The  analogy  of  the  shock  model 
in  risk  and  acturlal  analysis  has  been  given  by  BUhlmann  (1970,  Chapter  2). 

We  consider  Bayesian  nonparametric  procedures  for  the  shock  model  and 
assume  that  the  reader  is  familiar  with  the  results  in  Ferguson  (1973).  The 
prior  distribution  we  use  for  (X,F)  is  the  "independent  gamma-Dirichlet" 
prior.  By  this  we  mean  that  X  has  a  gamma  distribution  with  parameter 

Y  and  0  ,  denoted  by  U  a(dX)  ,  F  has  a  Dirichlet  distribution  with 

t»9 

parameter  a  ,  denoted  by  Pa(dF)  ,  and  X  and  F  are  independent.  We 

write  (X,F)  ~  U  a  *  P  and  find  that  the  posterior  distribution  of  (X,F) 
Y «“  « 

given  the  Poisson  process  up  to  time  T  and  Xj,  ....  X^^j  is  again  an 
Independent  gamma-Dirichlet  distribution,  namely,  ^y+{j(T)  0+T  *  P  * 

i 

where  the  summation  Is  from  one  to  N(T)  and  5x  denotes  the  probability 
measure  degenerate  at  x  . 


1 


1 

2 

The  aaymptotic  properties  of  the  posterior  distribution  is  considered 
and  it  is  shown  that  the  posterior  limiting  law  of  the  parameters  centered 
and  suitably  rescaled  is  Gaussian. 

A  parameter  of  primary  interest  is  the  survival  probability,  which 
is  the  probability  that  the  device  survives  beyond  time  t  .  We  calculate 
the  Bayes  estimate  for  this  parameter. 


2.  THE  PRIOR  AND  POSTERIOR  DISTRIBUTIONS 


Throughout  this  paper,  we  assume  that  the  damages  X^,  . ...  X^,  ... 
caused  by  the  shocks  are  independent  and  identically  distributed  random 
variables  having  distribution  F  and  that  the  shock  process 
Nt  *  (N(t)  ;  t  e  [ 0 , T] }  is  a  Poisson  process  with  parameter  X  which 
is  independent  of  the  damages.  The  observables  are  therefore  and 

h  "  (X1*  ••••  XN(T))  • 

We  formalize  the  above  assumptions  as  follows:  for  each  measurable 
partition  A^,  . . . ,  of  [0,T]  and  each  collection  ...  of 

Borel  sets  on  the  real  line,  the  joint  distribution  of  and  for 

given  (X ,F)  is  determined  by 


PX,F{NT(Aj) 


1. 


Xi£Ci 


i  -  1. 


(2.1) 


,  e-XTN(T) 


1 

n 

J-i 


ki 

«  3 (A,) 


V 


N(T) 

n  f(c  ) 

i-l  1 


N(T)  } 


where  N^(A^)  is  the  number  of  shocks  in  A^  ,  m  is  the  Lebesque 

measure  and  k^,  ...,  k^  are  nonnegative  integers  with 
l 


Ik-  N(T)  . 
j-1  3 


Using  the  idea  of  conjugate  priors,  we  let  X  ~  U  where  U 

y  >“  fi® 

2 

is  the  gamma  distribution  with  mean  y/6  and  variance  y/8  and  let 
F  where  Pa  is  the  Dirichlet  distribution  of  F  with  parameter 

a  where  a  is  the  finite  measure  on  the  reals.  With  the  additional 
assumption  that  X  and  F  are  Independent,  we  have  the  independent 
gaxmna-Dirichlet  distribution  Q  x  of  (X,F)  . 


We  need  the  following  notations  to  prove  the  main  theorem.  Let 
Qx  be  the  distribution  of  the  observations  given  X  ,  QF  the  distribution 


of  the  observations  given  F  and  Q  the  marginal  distribution  of  the 
observations.  The  characterizations  of  these  distributions  via  finite 
dimensional  sets  are  given  in  the  appendix.  We  are  ready  to  present  the 
following 


Theorem  2.1: 


Given 

N 

T 

and  X^.  ,  the  (posterior)  distribution  of 

(i) 

X 

is 

UY+N(T),0+T 

(ii) 

F 

is 

p  N(T)  and 

% 

' — 1 

o> 

jr 

i-l  i 

(iii) 

X 

and 

F  are  independent. 

Proof : 

Note  chat  according  to  (2.1),  given  (X,F)  ,  the  conditional  dis¬ 
tribution  of  Nt  and  Xj,  given  N(T)  and  Xp  is  independent  of  (X,F)  . 
Thus  {N(T),Xp}  is  a  set  of  sufficient  statistics  for  the  family  of  model 
distributions  F  :  E  R+  x®|  •  Thus,  we  can  treat  the  reduced 

problem  and  hence,  (2.1)  becomes 


(2.1’) 


Px  p{N(T)  »  k  ,  X±  e  Ct  ;  i  -  1 . N(T) } 


-XT/.-.k 

e  (XT) 

k! 


k 

n  f(c.)  . 
i-l 


To  prove  (i),  we  first  hold  F  fixed  (that  is  given  F) .  We  shall  show 
that  the  posterior  distribution  of  X  given  N(T)  and  X^  is 


(dX)  ,  which  is  independent  of  F  .  Ue  need  the  following 


Uy+N(T),0+T 


notations. 


k 

Let  A  -  {N(T)  -  k}  and  C  -  II  C  .  According  to  the  definition 

i-1 

of  conditional  distributions,  we  have  to  check  that  for  all  B  e  8  ,  for 
all  such  A's  and  for  all  such  C's  ,  the  following  equality  holds 


//\«UT>WB,Vdli(T>’d5>-,fPX,F<A  *  • 

AC  B 

k 

where  Q  is  given  by  (A.l).  Since  Q_(C)  -  II  F(C.)  is  cancelled 
F  F  ~  i-1  1 

out  in  both  sides  given  F  ,  we  need  only  check  that 


/  Uy+N(T),0+T(B)QF(dN(T)) 
A 


/ 


px,F(A)l\,e‘di>  • 


But  this  is  well  known  since  it  says  that  gamma  distributions  are  con¬ 
jugate  priors  when  sampling  from  Poisson  random  variables.  Thus  (i) 
holds  and  given  F  the  posterior  distribution  of  X  is  independent 
of  F  . 

Note  that  (ii)  can  be  proved  similarly.  We  first  hold  X  fixed, 

and  check  that  P  (dF)  satisfies  the  definition  of  the  posterior 

<*+•  V  6 
i-1  Xi 

distribution  of  F  given  N^  ,  and  X  .  We  need  to  check  if  the 
following  equality  holds  for  all  D  e  A  ,  all  A's  and  all  C's  , 
where  is  given  by  (A.l) 


/  /  P  N/T)  (D)Qx(dN(T).dx)  -/px>f(A  x  C)Pa(dF)  . 


A  -T  a +  l  & 

i-1  xi 


Since  QX(A) 
check  if 


-AT.,,_,k 
e  (XT) 

k! 


is  cancelled  out  in  both  sides,  we  need  only  to 


/P  k  (D)Qx(dx) 

5  a+  7  5 
i-1  Xi 


PX,F(?)Pa(dF) 


But  then  this  is  true  by  Ferguson  (1973,  Theorem  1  of  Section  3). 

We  have  also  shown  that  given  A  ,  the  posterior  distribution  of  F 
is  independent  of  A  ,  thus  (iii)  is  true. | | 


3. _ LIMITING  DISTRIBUTION 


The  consistency  of  the  posterior  distributions  is  well  known.  A 
general  proof  of  this  phenomenon  is  given  by  Doob  (1949)  who  uses  martingale 
methods  to  obtain  the  result  that  the  posterior  distribution  converges  to 
the  distribution  degenerate  at  the  true  parameter  under  which  the  sample 
is  actually  drawn,  for  almost  all  true  parameters.  The  delicate  problem 
of  the  classification  of  the  true  parameters  under  which  consistency  holds 
is  treated  by  many  authors  in  various  cases.  In  particular,  if  the 
parameters  are  discrete  distribution  functions  then  the  convergence  holds 
for  all  true  parameters  lying  in  the  support  of  the  so-called  tail  free 
prior  probabilities,  see  Freedman  (1963)  and  Fabius  (1964).  Berk  (1970) 
investigates  this  problem  in  the  dominated  case. 

In  this  section,  we  would  rather  search  for  the  limiting  posterior 
distribution  for  the  parameters.  Note  that  in  the  parametric  case, 
this  problem  has  been  considered  by  many  authors.  See  LeCam  (1963), 
or  more  recently,  Walker  (1969)  and  Dawid  (1970).  Our  problem  is  non- 
parametric  in  nature  and  our  results  are  based  on  a  theorem  of  Lo 
(1978a). 

We  denote  the  "true"  values  of  the  intensity  A  and  the  damage 
distribution  by  A  and  F_  respectively.  Thus,  P.  „  is  the  true 

0  o  aq,f0 

probability  distribution  of  the  observations.  We  denote  the  "posterior 
random  variable"  A  and  the  "posterior  Dirichlet  process"  F  given 
observations  up  to  time  T  by  A^  and  respectively.  Thus, 


8 


aT,FT)  ~ 


Uy+n(t),9+t  *  p  N(T) 


a+  7  <5 

U  % 


The  following  proposition  shows  the  consistency  of  the  posterior 
distributions. 

Proposition  3.1; 

The  posterior  distribution  of  X  converges  to  the  distribution 
degenerate  at  X^  and  the  posterior  distribution  of  F  converges  to 
the  distribution  degenerate  at  F^  ,  as  T  goes  to  infinity, 

‘•8-  [PV*o]  • 


Proof : 


Note  that  EUt  |  -  V/ffi  and  eJ|xt  -  I  Nj.XjJ 

.  Because  -*■  X  a.s.  IP.  _  1  via  some  renewal  arguments, 

T)2  T  0  L  X0'F0j 


y  +  N(T)  <  Because  -*■  X  a.s.  [p.  _  1  via  some  renewal  arguments, 

(0  +  T)2  T  0  L  0’F0j 

see  Chung  (1974,  p.  134).  Thus  the  posterior  mean  of  X  goes  to  X^ 
and  the  posterior  variance  of  X  goes  to  zero,  thus  the  first  assertion 
follows.  To  show  the  second  one,  let  B^  ;  j  *  1,  ...»  K  be  a  measurable 
partition  of  R  ,  we  shall  show  that  a.s.  P.  _  I  . 

L  vFoj 


E{Fj(Bj )  |  Nt,Xt>  -  FQ(B  )  ;  j  -  1 . K 


(2)  E<|Ft(Bj)  - 


N(T)  \ 

°+  l  6x  (M 

-  i-1  iJ  J  ' 

a(R)  +  N(T)  j 


/  N(T)  \ 

I  a  +  l  5X  (B  )1 
[F  (B  )  -L  i=L 

I  TV®J l'  a(R)  +  N(T)  j 


J  .  *  -  1.  ....  K 


Nt,Xt>^0  ; 


9 


Then  the  result  will  follow  via  Kallenberg  (1976,  Theorem  4,3),  Now, 
(1)  holds  because  as  T  goes  to  Infinity, 


E{FT(Bj) 


NT’V 


a^) 


N(T) 

+  l  *x  (B  ) 
i-1  xj  J 


a(R)  +  N(T) 


F0<»J> 


j  -  1 . K 


(2)  holds  because  for  j  l 


B<IVV  - 


N(T) 

a(V  +  Ji  5jti<Bj)  II  _  ”'*’1 

o(R)  +  N(T)  /I  W  a(R)  +  N(T) 


N(T) 

a(B.)  +  l  6  (B.)l 


i-1  Xi '  1 


E{FT(Bj)FT(BJl) 


w 


“<Bj)  + 


N(T) 

l  * 

i-1 


a(R)  +  N(T) 


E{W 


|  nt,xt} 


N(T) 

a(B  )  +  l  6  (B  ) 
i-1  i 

E(FT(Bj)  |  Rj.Xj} 

u(R)  +  N(T) 

N(T) 

“(V  +  1  «_  (B.) 

1  i-1  Xi  1 

N(T) 

a(B  >  +  I  «  (B.) 

J  i-1  i  3 

a(R)  +  N(T) 

a(R)  +  N(T) 

/ 

( °(b. )  +  t  &  (b,: 

\  J  i-1  xi  J 

’)(a<v + r  v^) 

(a(R)  +  N(T))(o(R)  +  N(T)  +  1) 

N(T) 

°(B. )  +  t  6  (B  ) 

3  i-1  i  3 

N(T) 

«(B.)  +  l  «  (B.) 

±ml  1  ■  n 

a(R)  +  N(T) 

a(R)  +  N(T)  U 

j  *  *  ,  j  ,  l 

*  If  « • • j  K  • 

The  case  for  j  -  A  can  be  proved  similarly.  |  | 


In  the  following,  ve  shall  prove  that  centered  and  properly  rescaled. 


Che  posterior  joinC  distribution  of  the  parameters  (X,F) 
weakly  to  a  joint  Gaussian  distribution,  a.s. 
be  distributions  on  [0,1]  .  Define 


[Vo] 


converges 
We  let  F 


N(T) 

o(t)  +  L  <$„  (t) 


Jtptt) 


»fo<f)  <FT(t)  - 


i-1  Xi 


a(l)  +  N(T) 


;  t  e  [0,1] 


and 


y  a  /f  I  x  -  *  n(t)  I 

*  tT  ■*'  )■ 

Let  Wq  ,  t  e  [0,1]  be  the  Brownian  bridge  subject  to  a  change  of  time 
by  Fq  ,  i.e.. 


EWQ(t)  •  0  ,  t  c  [0,1] 

E{WQ(t)  -  WQ(s)}2  -  {FQ(t)  -  FQ(s) Hi  -  FQ(t)  +  Fq(s) }  ,  s  <  C  . 

Let  Y  be  a  centered  normal  random  variable  with  variance  .  Then 

Theorem  3.1: 

If  F  is  continuous,  a.s.  P,  _  I 

0  L  VF0j 

(Yt,Xt(-))  (Y,W0(-))  when  T  -  - 

where  Y  and  Wg(*)  are  Independent. 


Proof; 


Because  Y^  and  X^,( • )  are  independent,  we  only  need  to  show 
(1)  YT— -Y  and  (2)  X^-)  -^WQ(*)  a.s.  ^  J  .  To  show  (1) , 

let  Q.  *  \N  :  -►  ;  .  This  set  is  of  P  probability  1;  see 

1  (  T  °!  VF0 

previous  arguments.  Our  assertion  will  follow  if  we  show  that  deleting 


a  P  -  null  set  and  for  all  N  c  0,  ,  and  all  {T  }  ,  T  ■*■<*> 

A0,F0  1  n  n 


we  have  Y_ 


•N(0,An)  as  n  -*■  »  .  See  Billingsley  (1968,  p.  16). 


Thus  we  pick  any  N  e  8^  ,  hold  N  fixed  and  then  let  T^  be  any  sequence 


of  positive  reals  that  goes  to  infinity.  We  shall  show  that  YT 

N(T  )  n 

Note  that  Y  -  J  (X  -  EX.)  where  X,  ;  j  -  0,1,  .... 
n  i-0  J 


N(0,A0) 


N(Tn) 


are  independent  gamma  random  variables  with  distributions  U 


V9+Tn 


respectively  with  y  •  y  ,  y  »  1  ;  1  <  i  <  N(T  )  .  Thus,  the  central 

u  x  —  —  n 


limit  theorem  for  triangular  array  applies  and  (1)  follows. 

(  N(T) 


We  show  (2)  similarly.  Let  8- 


.  ;n(t) 


«  (*) 


“  X 

i-1  xi 


•F  (•) 

0V  i 


This  set  is  of  P.  probability  1,  see  Proposition  3.1.  Deleting  a 

VF0 

P,  -  -  null  set,  fix  (N,X1 ,  ...,  X  ,  ...)  e  8«  and  let  T  be  a  sequence 
VF0  1  n  Z  n 

of  positive  reals  that  goes  to  infinity.  We  need  to  show  XT  — ' -  —  Wq ( » )  , 

n 

see  Billingsley  (1968,  p.  16).  But  then  this  is  proved  in  Example  1,  Part  II 
of  Lo  (1978).  Thus  a.s.  jp^J  ,  (X^O))  (Y,WQ,  (•))  ,  where 

Y  and  Wq(*)  are  independent. | | 

The  following  corollary  provides  an  asymptotic  Bayesian  simultaneous 
confidence  band  for  the  continuous  true  distribution  functions. 


Corollary  3.1: 


Under  Che  assumptions  of  Proposition  3.1,  we  have  a.s. 


lim  P  |  sup  jx_(t) 1  >  l!  -  2  l  (-1) 
T-*~  (0<t<l  )  j-1 


J+Ie-2^X2  ,  X  >  0 


I  \  -2X2 

lim  P  \  sup  XT(t)>X>»e  ,  X  >  0  . 

T-*~  ( 0<  t±l  ‘ 


Proof : 

These  are  well  known  consequences  of  Theorem  3.1.  For  detail 
arguments,  see  Billingsley  (1968,  p.  142). || 
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APPLICATIONS 


Lee  f  be  a  real  valued  integrable  (or  positive)  function  of  (A,F)  , 
then  the  Bayes  rule  given  the  observations  with  respect  to  a  quadratic 
type  loss  function  is  E{ f (A , F)  |  N^X^}  .  For  different  standard 
functions  of  F  ,  the  Bayes  rules  have  been  computed  in  Ferguson  (1973). 


Our  concern  here  is  the  survival  probability  H(t)  ,  H(t)  •  [  P, 


-At 


k-0 


k! 


where  is  the  probability  that  the  device  survives  k  shocks  for  the 

period  [0,t]  .  The  is  a  deterministic  function  of  k  and  the 

capacity  or  threshold  of  the  device.  Thus  H(t)  is  the  probability  that 
the  device  survives  the  period  [0,t]  .  We  consider  two  cases. 


Case  (1) : 

In  the  cumulative  damage  model,  the  kth  shock  is  survived  by  the 
device  if  X^  +  •  •  •  +  X^  does  not  exceed  the  capacity  or  threshold  y 
of  the  device.  Note  that  in  this  case  P^  -  PtX^^  +  •••  +  £  y  |  N(t)  ■  k}  . 

This  model  has  been  considered  by  Cox  (1962)  and  Barlow  and  Proschan  (1975) 
among  others. 

Case  (2a) : 

k  k 

P  ■  n  P{X.  <_  y. )  ■  II  F(y  )  .  This  represents  the  case  where 
k  i-1  1  1  i-1  1 

there  is  a  threshold  which  changes  after  each  shock  occurs  and  y^f  ...» 
are  the  successive  threshold  levels.  This  model  is  discussed  by  Esary 
et  •!.,  (1973). 


I 
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Cage  (2b) ; 

_  ^ 

-  { F (y ) }  .  This  Is  a  special  case  of  (2a)  with 

■  •••  -  y^  ■  y  ,  a  fixed  threshold  level. 

To  compute  our  Bayes  rule  for  the  survival  probability,  let  the 
loss  function  be 


L(H,fi)  -  /  (H<t)  -  iuH^Kdt) 


where  W  Is  a  totally  finite  measure  on  [O,-)  .  Thus,  the  Bayes 
rule  with  respect  to  this  loss  function  will  be,  for  each  t  , 

H(t)  -  E{H(t)  |  fcpXf}  . 


Note  that  E{H(t)  |  N^X^}  can  be  calculated  using  the  following  properties 
for  the  posterior  distribution  of  the  parameters. 


Property  (la) : 


E{e‘XtXk  |  NftXj. j 


r(y  +  N(T)  +  k)  /  1  \k/  8  +  T  \Y4-N(T) 

r(y  +  N(T))  y0+T+t/\0+T  +  t/ 


Proof: 


This  Is  a  consequence  of  Theorem  2.1. 
Property  (lb) : 


E |f(k) (y)  |  Nt,Xt}- 


N(T)  k 

a  +  >  6  +  T  6 

1-1  X1  1-1  U1 

Kh)  i 

N(T)  j-1 

ct  +  )  6  +  l  6 

1-1  xi  1-1  U1 

(dUj) 

o(R+)  +  N(T)  +  k  j-1 

o(R+)  +  N(T)  +  j  - 

■  1 
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Proof : 

This  follows  from  Theorem  2.1  and  repeated  applications  of  Lemma  1 
in  Lo  (1978b). 


Property  (lc) 

n  F(y  )  | 

(i-1  1 

Proof ; 

This  is  from  Theorem  2.1  and  an  application  of  Lemma  1  in  Lo  (1978b). 


Property  (Id) : 


E{Fk(y)  |  Nt,Xt| 


r {a(R  )  +  N(T))r 


(  N(T)  ) 

v»(y)  +  ^  «x  (y)  +  k> 


r{a(R  )  +  N(T)  +  k}r<o(y)  + 


N(T)  ) 

l  «  (y)> 

i-i  xi  ) 


Proof ; 

This  is  the  ktk  moment  of  beta  random  variables.  Now,  it  is  easy 
to  see  that 


E{H(t)  |  Nt>Xt) 


-Xt.k 

e  X 


*  Jo  1  • 


by  the  monotone  convergence  theorem 
E{.-AtXk  |  •  {£  , 


since  P^  depends  on  F  only  and  the  fact  that  under  the  posterior 
distribution,  F  and  X  are  independent. 
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•  lc 

•  l  E{F(k)(y)  I  NT,XT^E{e‘UXk  |  Nj.X^  •  jjy  in  Case  (1) 

<a»  i  ^  j  Jq 

•  IQ  E;A  F(y^)  |  NT,XTj  E|e"XtXk  I  Nt,Xt)  •  £7  in  case  (2a) 

-  I  E|Fk(y)  |  NT,XT|E^e-XtXk  |  ^  •  ~  in  Case  (2b), 

and  the  conditional  expectations  are  given  by  Property  (la),  (lb),  (lc) 
and  (Id). 

4.1  The  Case  of  Randoa  Threshold 

Suppose  there  is  no  practical  way  to  inspect  the  device  to  determine 

its  threshold  y  or  y^,  ....  y^ .  The  threshold  must  be  regarded 

as  a  random  variable.  We  let  y  —  G(dy)  and  in  case  of  y^,  ...»  y^  , 
we  let  these  be  independent  and  identically  distributed  as  G  .  The 
distribution  function  G  is  a  given  threshold  distribution  for  the  device. 

Of  course,  it  might  be  chat  G^dy^)  ;  i  -  1 . k  and  y^s 

are  independent.  But  this  case  can  be  treated  similarly  and  we  omit  the 
details  here.  It  is  not  difficult  to  see  that  Case  (1),  Case  (2a)  and 
Case  (2b)  become 


J  F(k)(y)G(dy)  . 

V 


and 
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Case  (2b'): 


?k  -  J  Fk(y)G(dy)  . 
0 


Now  it  is  easy  to  see,  using  Fubini's  theorem,  that  in  this  model, 

e|J  F(k)  (y)G(dy)  |  H^j  -  Je^F00^)  \  N^X^GUy) 

/»  00 

jf-fi  F(V  G<%>  l*rJrj 


00  00 


F(yi>  1  v?r|  j,  ^ 


eJJ  Fk(y)G(dy)  |  -  J  E^y)  | 

(o  o 


G(dy)  , 


and  the  Integrands  at  the  right  sides  are  available  in  Property  (lb),  (lc) 
and  (Id) . 
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APPENDIX 


Proposition  1: 


Qf{Nt(Aj)  -  kj  ,  J  -  1,  ....  1  ;  X±  e  ,  i  -  1,  ....  N(T) } 


(A.  1) 


Hi  ±  N(T»  /  J.  \Y(  1  f(T)  .  $  °  J  }  .  Nf }  F(c  ) 

T(y)  \  8  +  T/  \9  +  T/  A  kj  A  F(Ci}  • 


j-1  j 


i-l 


Qx{NT(Aj)  "  ^  *  1* 


A  ;  X  e  Cj.  ,  i  -  1,  ....  N(T) } 


(A. 2) 


«-XT(AT)N«>  n 


A  m  J(A. ) 


k  ! 

j-1  j 


/ 

N(T) 

n  c 

i-l  3 


NCT) 

n 

i-i 


i-l  1 

a  +  1  6X 

.  j-1  j- 

(dxt) 

o(R)  +  i  -  1 


T  '  j  '  j 


'i  & 


(A.  3) 


r(y  +  N(T))  /  8  \yl  1  \N(T)  J  m  (V  f  L  1=1  XjJ  1 

r(Y)  \8+T/\e+T/  A  k,!  J  A  a(R)+i-l 

J  A  J  MfT'S  1  A 


1 

.  N(T) } 

'  i-l 

a  +  I  <5 
.  j-1  Xj_ 

(dx£) 

a(R)  +  i 

-  1 

N(T) 

E  C 
i-l  3 


Proof ; 

(A. 1)  is  obtained  by  integrating  (2.1)  with  respect  to  U  .  To 

Y*° 

show  (A. 2)  and  (A. 3),  we  first  integrate  (2.1)  with  respect  to  P^  and 

U  x  P  respectively  and  then  an  application  of  Lemma  1  in  Lo  (1978b) 
Y*  “  o 

concludes  the  proof. || 


